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Abstract 

We discuss the resummation approach in QCD Analytic Perturbation Theory (APT). We start with a simple exam- 
ple of asymptotic power series for a zero-dimensional analog of the scalar g ip 4 model. Then we give a short historic 
preamble of APT and show that renormgroup improvement of the QCD perturbation theory dictates to use the Frac- 
tional APT (FAPT). After that we discuss the (F)APT resummation of nonpower series and provide the one-, two-, and 
three-loop resummation recipes. We show the results of applications of these recipes to the estimation of the Adler 
function D(Q 2 ) in the N/ = 4 region of Q 2 and of the Higgs-boson-decay width F H ^(m^) for M H = 100-180 GeV 2 . 

Keywords: Renormalization group, QCD, Analytic Perturbation Theory, Nonpower Series Resummation, Adler 
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1. Simple example of asymptotic power series 

In spite of many examples of successful applications 
of perturbative approach in quantum field theory, a per- 
turbative power expansion for a quantum amplitude usu- 
ally is not convergent. Instead, a typical power series 
appears to be asymptotic. To refresh the reader knowl- 
edge on this subject we consider a simple example — 
the so-called "0-dimensional" analog of the scalar field 
theory gip 4 : 



/-» CO 

1(g) = e-x-^dx. 

%J—oo 

It can be expanded in a power series [Q| 



(la) 



(lb) 



k=0 



with factorially growing coefficients: 4 = F(2k + 
l/2)/Y(k + 1) 2 k k\ fork » 1. Meanwhile, 1(g) 
can be expressed via special MacDonald function 



Kg) 



1 



2 V£ 



Ki/4 I — 



(lc) 



with known analytic properties in the complex g plane: 
It is a four-sheeted function analytical in the whole com- 
plex plane besides the cut from the origin g = along 
the whole negative semiaxis. It has an essential singu- 
larity e~^ 8s at the origin and in its vicinity on the first 
Riemann sheet it can be written down in the Cauchy in- 
tegral form: 



i(g) = VSF- 



V2~tt 



n Jo 



dye 



-l/4y 



jig + r) 



(id) 



•Responsible for Sect. l4.3l andl5l 



Due to this singular behavior near the origin the power 
Taylor series (ITbb has no convergence domain for real 
positive g values. This behavior is in one-to-one corre- 
spondence with the factorial growth of power expansion 
coefficients. The same factorial growth of expansion co- 
efficients has been proved for the <p 4 scalar and a few 
other QFT models M- 

Power series with factorially growing coefficients be- 
longs to the class of Asymptotic Seria (AS) — their 
properties were investigated by Henry Poincare at the 
end of the XIX century. In short, he concluded that the 
truncated AS can be used for obtaining the quantitative 
information on expanded function. To be more con- 
crete, the error of approximating F(g) by first K terms 
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Table 1: The last detained ((-g) Ik) and the first dismissed ((-g) 1k+\) contributions to the power series 1 1 bt in comparison with the exact 
value 11c) and the approximate result Ik(s), with K being the truncation number and A^Iig) — the error of approximation. 



g 


K 






I K (g) 


Kg) 


A*/(s) 


0.07 


7 


-0.04(2%) 


+0.07(4.4%) 


1.674 


1.698 


1.4% 


0.07 


9 


-0.17(10%) 


+0.42(25%) 


1.582 


1.698 


7% 


0.15 


2 


+0.13(8%) 


-0.16(10%) 


1.704 


1.639 


4% 


0.15 


4 


+0.30(18%) 


-0.72(44%) 


1.838 


1.639 


12% 



of expansion, F(g) -> F K (g) = Z k <Kfk(g), is equal to 
the last accounted term //Kg)- This observation can be 
used to obtain a lower limit of possible accuracy for the 
given g value: one should find the number K with the 
minimal value of ficig) — then its absolute value just 
set the limit of accuracy. 

To make this statement more clear, we take a power 
AS, fk(g) = fk g k , with factorial growth of expansion 
coefficients, fk ~ k\. Then we see that terms fk(g) stop 
to decay at k ~ K ~ l/g, hence for k > K + 1 trun- 
cation error starts to grow. The explicit calculations for 
the function Kg) ( flat with AS (fTBl . presented in TableQ] 
(see also in [3J), shows that the optimal values of trun- 
cation number K — K* l/(2g) (in dark green) indeed 
provides the best possible accuracy (dark green in the 
last column): Attempting to account a couple of extra 
terms (the second and the forth lines) results in drastic 
rise of the error! Thus, one has K,(g - 0.07) = 7 and 
K,(g - 0.15) = 2 . It is not possible at all to get the 1% 
accuracy for g = 0.15 , instead the best accuracy here is 
4%. In this situation to obtain more accurate result we 
need to use Eq. ([Tel l which has the right analytic proper- 
ties in g. 

We conclude from this example that to produce reli- 
able predictions for a physical amplitude in QCD one 
needs a tool for resumming the original perturbative 
power series into the quasi-nonperturbative amplitude 
which should have prescribed analytical properties, dic- 
tated by the first principles of QCD, like causality and 
renormalizability. It is intriguing that the APT approach 
in QCD allows one to obtain such a tool. 



2. Analytic Perturbation Theory in QCD 

In the standard QCD PT we know the Renormaliza- 
tion Group (RG) equation in the /-loop approximation 



da { i)[L] 
dh 



1 +Y J c k a\ l) \L] 



k>] 



for the effective coupling a s (Q 2 ) = a(r,[L]//3f with L = 
ln(g 2 /A 2 ), p f = b (N f )/(4n) = (11 - 2N f /3)/(4x)R 
Then its one-loop solution generates Landau pole sin- 
gularity, fl(i)[L] = l/L. As a consequence, the pertur- 
bative power series for the Adler function, D(Q 2 ) = 
do + 2jt>i dka k [L], being reasonably good in the deep 
UV region of Q 2 , at Q 2 = A 2 has the unphysical pole 
— in marked contrast to the analyticity property of the 
Adler function. Indeed, it should have a cut along the 
negative axis of Q 2 , but can not have any poles in Eu- 
clidean domain, Q 2 > 0. APT was suggested as a res- 
olution of this contradiction: From the very beginning 
we demand that improved coupling, as well as all its in- 
teger powers, should have the right analytic properties, 
i. e. be represented as dispersive integrals, see ©. 

Strictly speaking the QCD Analytic Perturbation 
Theory (APT) was initiated by N. N. Bogoliubov et 
al. paper of 1959 J3], where ghost-free effective cou- 
pling for QED has been constructed. Then in 1982 
Radyushkin ^ and Krasnikov and Pivovarov |0] us- 
ing the same dispersion technique suggested regular (for 
s > A 2 ) QCD running coupling in Minkowskian region, 
the well-known n~ l arctan(7r/L). After that in 1995 
Jones and Solovtsov using variational approach 10] 
constructed the effective couplings in Euclidean and 
Minkowski domains which appears to be finite for all 
Q 2 and s and satisfy analyticity integral conditions (0. 
Just in the same time Shirkov and Solovtsov l8|], us- 
ing the dispersion approach of discovered ghost- 
free coupling ^7li[L], Eq. (l5ab . in Euclidean region and 
ghost-free coupling 2ti [L], Eq. (l5bl . in Minkowskian re- 
gion, which satisfy analyticity integral conditions 

^(e^DtSliKe 2 )^ 2 f \i, 2 dcr; (3a) 

Jo 



%ds)=R[Jli](s) 



r 

2m J_ s 



(a- + Q 2 )- 
i+,£ ^i(cr) 



dev. (3b) 



(2) 



We use notations f(Q ) and f[L] in order to specify the argu- 
ments we mean — squared momentum Q 2 or its logarithm L = 
ln(Q 2 /A 2 ), that is f[L] = /(A 2 • e L ) and the QCD scale parameter 
A is usually referred to Nt = 3 region. 
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The last coupling coincides with the Radyushkin one 
for s > A 2 . Due to the absence of singularities in these 
couplings, Shirkov and Solovtsov suggested to use them 
for all Q 2 and s (for recent applications — see in |@]). 

Shirkov-Solovtsov approach, now termed APT, ap- 
pears to be very powerful: in Euclidean domain, -q 2 = 
Q 2 , L = \nQ 2 /A 2 , it generates the following set of 
images for the effective coupling and its n-th powers, 
{yt„[L]}„ eN , whereas in Minkowskian domain, q 2 = s, 
L s = In s/ A 2 , it generates another set, {2l„|X s ]} n6 N- APT 
is based on the RG and causality that guaranties stan- 
dard perturbative UV asymptotics and spectral proper- 
ties. Power series of the standard PT £ m d m a"'[L] trans- 
forms into non-power series £„, d m M m [L] in APT. 

By the analytization in APT for an observable f(Q 2 ) 
we mean the dispersive "Kallen-Lehman" representa- 
tion 



Pficr) 



cr+Q 2 - ie 



■ dcr 



(4) 



with pf(cr) = n 1 Im [/(— cr)]. Then in the one-loop 
approximation p\ (cr) = 1/ ^L 2 + n 2 and 



r pi(q-) 

Jo cr + Q 2 

f 



dcr - — - 



L e L -V 



Pi(<r) , 1 



dcr = — arccos ■ 



» ^7T 2 + L 2 



(5a) 



(5b) 



whereas analytic images of the higher powers (« > 
2, n G N) are: 



rtn\L]\ 1 / d\"- L (&i'[L] 



%\?[L S ]} (ft- 1)1 \ dLj \^ [Ls] 



n-l i rnWr 



(6) 



3. Fractional APT in QCD 

At first glance, the APT is a complete theory pro- 
viding tools to produce an analytic answer for any 
perturbative series in QCD. But in 2001 Karanikas 
and Stefanis II 1 011 suggested the principle of analyti- 
zation "as a whole" in the Q 2 plane for hadronic ob- 
servables, calculated perturbatively. More precisely, 
they proposed the analytization recipe for terms like 
j Q dxf Q dy a s (Q 2 xyj f(x)f(y), which can be treated as 
an effective account for the logarithmic terms in the 
next-to-leading-order approximation of the perturbative 
QCD. This actually generalizes the analytic approach 
suggested in 111 111 . Indeed, in the standard QCD PT one 
has also: 

(i) the factorization QCD procedure that gives rise to the 
appearance of logarithmic factors of the type: a v [L] L; 

(ii) the RG evolution that generates evolution factors of 



the type: B(Q 2 ) = [Z(Q 2 ) /Z(p 2 )] B{p 2 ), which reduce 
in the one-loop approximation to Z(Q 2 ) ~ a v [L] with 
v - 7o / (2£>o) being a fractional number. 
All that means that in order to generalize APT in the 
"analytization as a whole" direction one needs to con- 
struct analytic images of new functions: a" ', a v L m , 

This task has been performed in the frames of the so- 
called FAPT, suggested in IU2U13I1 . Now we briefly de- 
scribe this approach. 

In the one-loop approximation using recursive rela- 
tion © we can obtain explicit expressions for both cou- 
plings: 

F(e- L , 1 - v) 



= 77 



sin 



(v- 



T(y) 

1) arccos ( , = ) 

\ Vjr 2 +L 2 / 



n(v- 1)(tt 2 + L 2 ) 



2 \(v-l)/2 



(7a) 



(7b) 



Here F(z, v) is reduced Lerch transcendental function, 
which is an analytic function in v. The obtained func- 



tions, ^^[L] and 2ly'[L], have very interesting prop- 
erties, which we discussed extensively in our previous 
papers Il2l4l5ll . Note here, that in the one-loop approx- 
imation to find analytic images of aTjJ-L] • L m is very 

easy: they are just ^l m [L] and Sl^JL]. 

Constructing FAPT in the higher-(/)-loop approxima- 
tions is a more complicated task. Here we represent the 
original /-loop running coupling in the following form: 

JfmVLA 



AV) 



fl W [L a - in] = 



(8) 



with R(f) [L] and <p(Q [L] being the known functions. Then 
the spectral densities of the v-power of the coupling is 



(0 _ 1 sin[v f(i)[Lg-]] 
n (R (l) [La-]) 

but spectral densities of a v [L] ■ LP 
plicated: 



(9) 



start to be more com- 



P ( £ [L] = i'v , n sin - ™<p m [L]] . 

i " J " nR ([) [L] 



(10) 



Construction of FAPT with fixed number of quark fla- 
vors, Nf, is a two-step procedure: we start with the per- 
turbative result |a(C2 2 )] , generate the spectral density 
Pv(cr) using Eq. ©, and then obtain analytic couplings 
3i v [L] and 2l v [L] via Eqs. (0. Here Nf is fixed and fac- 
torized out. We can proceed in the same manner for 
A^-dependent quantities: [c^Q 2 ; iV/)J => p v (o~,Nf) = 

p v [La-,N f ] = pM)IB} => &y[L;Nf] and %[L;N f ] — 
here Nf is fixed, but not factorized out. 
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Then d„ - ((^'^p^ an d we need to resum the series 
nV v [L; r; T\ = J] r '"' ^[X] , (16a) 



Global version of FAPT 11411 . which takes into ac- 
count heavy-quark thresholds, is constructed along the 
same lines but starting from global perturbative cou- 
pling [a!j" >b ((2 2 )l , being a continuous function of Q 2 due 
to choosing different values of QCD scales A/, corre- 
sponding to different values of Nf. We illustrate here 
the case of only one heavy-quark threshold at s = rrv\, 
corresponding to the transition Nf = 3 — > Nf = 4. Then 
we obtain the discontinuous spectral density 

pfV) = 8(L cr <L 4 )p n [La-; 3] 

+ e(L 4 <L cr )p n [L . + A 4 ;4] , (11) 

with La- = ln^cr/Aj), Lj s ln^/n^/Aj) and Af = 

ln^Aj/A^J for / = 4, which is expressed in terms 
of fixed-flavor spectral densities with 3 and 4 flavors, 
p n [L; 3] and p„ [L + A 4 ; 4] . However it generates the con- 
tinuous Minkowskian coupling 

2l-* b m = 6(L<L 4 ) (ajL; 3] + A 43 2t v ) 

+ G(L 4 <L) 2t„[L + /U; 4] . (12a) 

with A 43 2t v = 2i v [L 4 + /1 4 ; 4] - 2l v [L 4 ; 3] and the analytic 
Euclidean coupling Sly 0b [L] 

pjglob 



serin = si v [L + A 4 -A] 



p 



p v [L <T ;3]-p v [La- + A 4 ;4] 



1 + e 



dLa- (12b) 



(for more detail see in 11411 ) 



4. Resummation in (F)APT 

Before starting with the definite-loop approximation 
we introduce here, following [16], the generating func- 
tion P{t) for perturbative coefficients c4 which allows us 
then to resum any non-power series (v = corresponds 
to the APT case) of the type 

S V [L\ T] = d Q T V [L) +d l J] d„ T n+V [L] , (13) 

n>\ 

where T\L\ denotes one of the analytic quantities 
m (,) [L\, 21® [L], or p (/) [L], and d„ = djdi. We suppose 
that 

Jr»oo /-*co 
P{t)t"- l dt with I P(t)dt=l. (14) 
Jo 

To shorten our formulas we use the abbreviated notation 
F[L,t]P(t)dt. (15) 





related to the original one in a simple way 

S V [U T] = d T V [L] + d x «nV v [L; t; T])) P(t) . (16b) 

4.1. One-loop FAPT 

In the one-loop approximation we have the following 
recurrence relation CF[L\ = dT\L\jdL): 

1 



' Jn+vin — ^Fn+l+vlL] 



(17) 



This property of couplings and spectral densities allows 
us to resum the series ( I16bl ): 

S V [L, n = d Q T V [L] + <?, ((T l+V [L - t])} PM , (18a) 

where now the generating function P v depends on v, 

dx 



Py(t) 



t>v(x) 



1-x 



(18b) 
6(x), and 



Here O v (jc) = vx y so that lim^o^v 
therefore lim v ^ P v (0 = P(0- 

4.2. Two-loop FAPT 

In the two-loop approximation we have more compli- 
cated recurrence relation 

1 



■ r„ + AL] = T n+ i +v [L] + c, T„ +2+ v[L] (19) 
n -1- v 

with c\ being the corresponding coefficient in Eq. ©. 
In order to resum the series 'WvlL; t\ T\ we need to in- 
troduce the "two-loop evolution" time 



T 2 (f) = t - c\ In 



1 + — 



; t 2 (0 = 



1 



(20) 



l+d/t 

We obtained in jl7] the following resummation recipe 

W t; T] = T v+ i [L] + A(v) ci t 2 (f) r 2 [L f , ] 
1 

-t 2 (0 J z v [f 7"v + i [L«J - ^ r v+2 [L u ]] dz (21) 


with L,,- - L + T 2 (fz) - T 2 (f), L,^ = L — T 2 (f), and 
A(v) being a Kronecker delta symbol. Interesting to note 
here that it is possible to obtain an analogous, but more 
complicated recipe for the case when T is the analytic 
image of the two-loop evolution factor a v (l + cia) v ' , see 
in [17] for more detail. 
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Table 2: Coefficients d„ for the Adler-function series with Nf = 4. The numbers in the square brackets denote the lower and the upper limits of the 
INNA estimates. 





PT coefficients 




ch 


d A 


d s 


1 


pQCD results with N f = 4 [18, 19) 


1 1.52 


2.59 


21 A 




2 + 


Model (26) with c = 3.544, S = 1.3252 


1 1.53 


2.80 


30.9 


2088 


2 


Model <26) with c = 3.553, S = 1.3245 


1 1.52 


2.60 


27.3 


2025 


2" 


Model {26) with c = 3.568, <5 = 1.3238 


1 1.52 


2.39 


23.5 


1969 


3 


"INNA" prediction of [17.1 


1 1.44 


[3.5,9.6] 


[20.4, 48.1] 


[674, 2786] 



4.3. Three-loop FAPT 

We describe here the recently obtained results on re- 
summation in the three-loop FAPT. In this case the re- 
currence relation has three terms in the r.h.s. 
1 



« + v 



+ C l T„ + 2 + v[L]+C 2 Tn + 3 + v[L] (22) 

with co being the corresponding coefficient in Eq. (|2). 
We introduce the "three-loop evolution" time by 

c]-2c 2 \ tA 

t 3 (0 



t + 



" 2 ln 



arc tan 



1 ci+t 
1 + 1 



2 C2 + C\ t 



(23) 



dr 3 (t) 



= hit) = 



l 



dt " J "' 1 +ci/t + c 2 /t 2 ' 
Then our resummation recipe is 

nv v [L, t; r] = r v+l m + am c 2 t 3 (0 r 3 [L,, ] 



+«f v+2 [L]-T 3 (() z 



z 



+zf T V+2 [L,J + (v + 1) f T v+2 [L f , : ] 



& (24) 



with L,^ = L + r 3 (f z) - r 3 (f) and L,$ - L - r 3 (f)- 



5. Resummation for Adler function 

Here we consider the power series of the vector corre- 
lator Adler function (labeled by the symbol V) Hi, [HI 



D y [L] = 1 + 



n>\ v ' 



(25) 



Due to d\ = 1 coefficients d„ coincide with d„. We 
suggested lfl7ll the model for the generating function of 
the perturbative coefficients d„ (see 1st row in Table [2) 

6e-' /cS -(t/c) e-' /c 



Py(t) 



(26a) 



which provides the following Lipatov-like coefficients 

5" +l - n 



d: = d 



n-\ 



T(n). 



(26b) 



6 2 - 1 

Our prediction c/J = 27.1, obtained with this gen- 
erating function by fitting the two known coeffi- 
cients d 2 and dj and using the model (126) . is in 
a good agreement with the value 27.4, calculated 
in Ref. Jlil Note that fitting procedure, tak- 

ing into account the fourth-order coefficient d\, pro- 
duces the readjustment of the model parameters in 
d26) to the new values {c = 3.5548, 5 = 1.32448} -> 
{c = 3.5526, 6 = 1.32453). The corresponding values 
of coefficients d^ are shown in the third row labelled by 
2 in Table |2] 

In order to understand how important are the exact 
values of the higher-order coefficients d„, we employed 
our model ( 126) with two different sets of parameters 
c and 6, shown in rows labelled by 2 + and 2~ in Ta- 
ble [2] One set, 2 + , roughly speaking, enhances the ex- 
act values of the coefficients c/ 3 and c/4 by approximately 
+8% and +13%, correspondingly, while the other one, 
2~, — reduces them in the same proportion. All coef- 
ficients of these models are inside the range of uncer- 
tainties determined in Il7ll using the Improved Naive 
Non-Abelinization (INNA). Moreover, the difference 
between the analytic sums of the two models in the re- 
gion corresponding to Nf =4 is indeed very small, 
reaching just a mere +0.05%. This gives an evident sup- 
port for our model evaluation. 

Now we are ready to estimate the relative errors, 
A^[L], of the APT series^] truncation at the Mh term: 



D. 



n>l 

Dlm-D V N [L] 



K[L] 



Dim 



(27a) 
(27b) 



2 Note that power series 1251 has v = — for this reason we use 
here the APT approach. 
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Here is the resummed APT result in the cor- 

responding /-loop approximation, see Eqs. ( 1 1 8ab . ( f2Tb . 
and (l24l) with substitution v — > 0. In Fig. [T] we show 
these relative errors for N - 1, 2, 3, for the one- and 
two-loop cases (calculations for the three-loop case is 
not yet finished). The main result is in some sense sur- 
prising: The best order of truncation of the FAPT series 
in the region Q 2 = 2 - 20 GeV 2 is reached by employ- 
ing the N 2 LO approximation, i.e., by keeping just the 
c/2-term. 



0.014f\ 
. 012 
0.01 
.008 
0.006 
. 004 







A^ 


Q 2 [GeV 2 ] - 



2.5 5 7.5 10 12.5 15 17.5 20 

Figure 1: The relative errors AjJ(g 2 ) evaluated for different values of 
N: N = 1 (short-dashed red line), N = 2 (solid blue line), and N = 3 
(dashed blue line) of the truncated APT given by Eq. (27b) , in com- 
parison with the exact result of the one- and two-loop resummation 
procedure represented by Eqs. 1 1 8i and )2 1 \ . 

We may also compare the numerical values for the 
resummed quantities, obtained in different loop approx- 
imations. We take for this comparison the following l- 



loop QCD scale parameters at Nf = 3 flavors: A 



201 MeV, A 



(1=2) 



f 

379 MeV, and A 



(i=i) 



(1=3) 



385 MeV, 



which have been determined from the condition that the 
APT prediction for the ratio R e + e -(s = m|) should coin- 
cide with the "experimental" value 1 .03904, determined 
in IU9I1 . We obtain the following values of the resummed 
Adler functions, shown in the table form for two values 
of Q 2 , namely 3 and 2 GeV 2 : 



Loop order / 


1=1 1=2 1 = 3 


£^ (0 (e 2 = 3GeV 2 ) 


1.1129 1.1131 1.1164 


Dl (,) (Q 2 = 2 GeV 2 ) 


1.1221 1.1223 1.1257 



Strictly speaking, our model for coefficients is valid 
only for Q 2 > m] = 2.44 GeV 2 , so that the first value 
Q 2 = 3 GeV 2 was selected to show the results in the 
legitimate Nj = 4 domain. The second value, Q 2 = 
2 GeV 2 , was selected for the comparison with the recent 
estimate in Ir20ll . where the value DZ,(Q 2 ) = 1.1217 has 
been obtained in the two-loop approximation using the 
so-called generalized Pade summation method. 



6. Resummation for H° -* bb Decay Width 

Here we analyze the Higgs boson decay to a bb pair. 
For its width we have 



r(H = -^M H R S (M 2 H ) 
4V2tt 



(28) 



with R S (M 2 H ) = m 2 b {M 2 H )R s {M 2 H ) and R s (s)is the fl-ratio 
for the scalar correlator, see for details in Il2l 12111 . In the 
one-loop FAPT this generates the following non-power 
expansion: 

K S [L] = 3 m 2 w U^[L] +dlY J '^ - (29) 

where fiit\) = 8.21 - 8.53 GeV is the RG-invariant of 
the one-loop nibifi 2 ) evolution m 2 b {Q 2 ) = m 2 ^a v °(Q 2 ) 
with vo = 2yo/£>o(5) = 1.04 and -yo is the quark-mass 
anomalous dimension (for a discussion — see in lfl7ll )PI 
We take for the generating function P{t) the model 
of JH with {c = 2.4, p = -0.52} 



c(l +P) 

It provides the following Lipatov-like coefficients 
n _iT{n + \)+f}T{n) 



d! = c n 



(30a) 



(30b) 



which are in a very good agreement with d*,n = 2, 3, 4, 
calculated in the QCD PT 12 ill , the corresponding val- 
ues of coefficients d* are shown in the third row labelled 
by 2 in Table [3] In order to estimate the importance 
of the higher-order coefficients d„ exact values, we pro- 
ceed along the same lines as in Sect. [5] We employ our 
model (f30b with two different sets of parameters c and 
yS, shown in rows labelled by 2 + and 2" in Table [3] One 
set, 2 + , enhances the exact values of the coefficients 
and d\ by approximately +13% and +20%, correspond- 
ingly, while the other one, 2~, — reduces them in the 
same proportion. Obtained in this way difference be- 
tween the analytic sums of the two models in the re- 
gion corresponding to M# = 80 - 170 GeV is indeed 
very small, not more than 0.5%. Note here also that 
the model prediction for c/5 is very close to the predic- 
tion obtained using the Principle of Minimal Sensitivity 
(PMS) JH, shown in the row with label 3 in Table |3] 



3 Different values of ffi ( i) is related with two different extractions 
of the RG-effective fc-quark mass, mt(nr b ), which we have taken 
from two independent analyses. One value originates from Ref. 1241] 



(mi,(mp=4.35: 
yielding mi,(mf) 



0.07 GeV), while the other was derived in Ref. (25 
= 4.19 ±0.05 GeV. 
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Table 3: Coefficients d„ for the Higgs-boson-decay width series with Nr = 5. 





PT coefficients 




d 3 




d s 


1 


pQCD results with N f = 4 [211 


1 7.42 


62.3 


620 




2 + 


Model OS) with c = 2.43, <S = -0.52 


1 7.85 


68.5 


752 


10120 


2 


Model OS} with c = 2.62, 6 = -0.50 


1 7.50 


61.1 


625 


7826 


2" 


Model 0S1> with c = 2.25,(5 = -0.51 


1 6.89 


52.0 


492 


5707 


3 


"PMS" predictions of [21-231 




64.8 


547 


7782 



After verification of our model quality we apply the 
FAPT resummation technique to estimate how good is 
FAPT in approximating the whole sum KslL) in the 
range L e [12.4, 13.5] which corresponds to the range 
Mh e [100, 170] GeV 2 with A^™ = 201 MeV. Here 
we need to use our resummation recipes ( 1 1 8ab and 
(l2lT) with substitution v -> v = 1.04. Note that in 
the two-loop approximation the one-loop evolution fac- 
tor a v transforms into a more complicated expression, 
a v (l + c\ a) v ' . This produces additional numerical com- 
plications, but qualitatively results are the same. For this 
reason we show explicitly only one-loop formulas. 

We analyze the accuracy of the truncated FAPT ex- 
pressions 



N 3s 



^m+di 2^ sew 



(31) 



and compare them with the resummed FAPT result 
"R S [L] in the corresponding /-loop approximatiorfl using 
relative errors A N [L] = l-'R s [L;N]/'R s [L]. We estimate 
these errors for N = 2,N = 3, and N - 4-in the analyzed 
range of L e [1 1, 13.8] and show that already ^ S [L; 2] 
gives accuracy of the order of 2.5%, whereas < R S [L;3] 
of the order of 1%. That means that there is no need to 
calculate further corrections: at the level of accuracy of 
1% it is quite enough to take into account only coeffi- 
cients up to ds. This conclusion is stable with respect to 
the variation of parameters of the model P s (t) and is in a 
complete agreement with Kataev-Kim conclusion 12711 . 

In Fig. |2] we show the results for the decay width 
r^^(Mtf) in the resummed two-loop FAPT, in the 
window of the Higgs mass allowed by the LEP and 
Tevatron experiments. Comparing this outcome with 
the one-loop result, shown as the upper strip in Fig. [2] 
reveals a 5% reduction of the two-loop estimate. This 
reduction consists of two parts: one part +7%) 
comes from the difference in the mass m, while the other 




100 



110 



120 



130 



140 



150 



160 



170 



Figure 2: The two-loop width . in the resummed FAPT as a 
function of the Higgs-boson mass Mfj . The mass is varied in the in- 
terval rh(2) = 8.22 ± 0.13 GeV according to the Penin-Steinhauser 
estimate 12411 . The upper strip shows the corresponding one-loop re- 
sult. 



-2%) is due to the difference in the values of Rs (Mh) 
in both approximations. 

7. Conclusions 

We conclude with the following. 

APT provides natural way to Minkowski region for 
coupling and related quantities with weak loop depen- 
dence and practical scheme independence. 

FAPT provides an effective tool to apply the APT ap- 
proach for the renormalization-group improved pertur- 
bative amplitudes. 

Both APT and FAPT produce finite resummed an- 
swers (now — up to the three-loop level) for perturba- 
tive quantities if we know the generating function P{f) 
for the PT coefficients. 

Using quite simple model generating function P(t) 
for the Adler function D(Q 2 ) we show that already at 
the N 2 LO an accuracy is of the order 0.1%, 
whereas for the Higgs boson decay H — > bb at the N 3 LO 
is of the order of: 1 % — due to the truncation error; 2% 
— due to the RG-invariant mass uncertainty. 



4 Here we show the results only for 1=1 and / = 2: Calculations 
for 1 = 3 are not yet finished. 
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